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Abstract 

A new constant of motion for helically symmetric equilibria in the vicinity 
of the magnetic axis is obtained in the framework of Vlasov theory. In view of 
this constant of motion the Vlasov theory is compared with drift kinetic and 
gyrokinetic theories near axis. It turns out that as in the case of axisymmetric 
equilibria [H. Tasso and G. N. Throumoulopoulos, Phys. Plasmas 18, 064507 
(2011)] the Vlasov current density thereon can differ appreciably from the 
drift kinetic and gyrokinetic current densities. This indicates some limitation 
on the implications of reduced kinetic theories, in particular as concerns the 
physics of energetic particles in the central region of magnetically confined 
plasmas. 
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I. Introduction 



Kinetic equilibria may provide broader and more precise information than 
multifluid or magnetogydrodymanic equilibria as those governed by the Grad- 
Shafranov equation. Since solving self consistently the kinetic equations is 
tough particularly in complicated geometries the majority of kinetic equilib- 
rium solutions are restricted to one dimensional configurations in plane ge- 
ometry, e.g. PQ-[E]- Of particular interest are equilibria with sheared toroidal 
and poloidal flows which play a role in the transition to improved confine- 
ment regimes in tokamaks and stellarators, though understanding the physics 
of this transition remains incomplete. Construction of kinetic equilibria is 
crucially related to the particle constants of motion which the distribution 
function depends on. In the framework of Maxwell- Vlasov theory only a 
couple of constants of motions are known for symmetric two-dimensional 
equilibria, i.e. the energy, H, and the momentum p x s conjugate to the ignor- 
able coordinate, £3, out of the four potential constants of motion. Therefore, 
for distribution functions of the form f(H,p x3 ) only macroscopic flows and 
currents along the direction associated with x 3 can be derived, e.g. toroidal 
flows for axisymmetric plasmas. The creation of poloidal flows requires addi- 
tional constant (s) of motion. This remains an open question despite the fact 
that in a previous study [9] we found locally in the vicinity of the magnetic 
axis the following new constant of motion: C = v z + Jin \v^\, where is the 
toroidal particle velocity, v z the velocity component parallel to the axis of 
symmetry and = I/re^ the toroidal magnetic field near axis (r,(f),z are 
cylindrical coordinates). An additional open question remains a potential 
extension of the proof of non existence of magnetohydrodynamic axisym- 
metric equilibria with purely poloidal flows [TO] to Vlasov equilibria. A third 
constant of motion was studied extensively in the Astrophysics literature 
[TT]-[21]. In particular for axisymmetric astrophysical systems distribution 
functions depending on an approximate third constant of motion result in 
velocity ellipsoids with unequal axes in agreement with observational data 
for our Galaxy O EH OS] . 

Owing to the contemporary and probably future limited computational 
efficiency for simulations in the framework of Vlasov theory approximate 
kinetic theories have been established in reduced phase spaces as the drift 
kinetic and gyrokinetic ones. In both theories the reduced phase space is 
five dimensional with three spatial components associated with the guiding 
center position, R (or the gyrocenter position in the framework of gyrokinetic 
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theory), and a velocity component parallel to the magnetic field, v»; also, 
the two components of the perpendicular particle velocity are approximated 
after a gyroangle averaging with the magnetic moment which is treated as 
an adiabatic invariant. A related underlying assumption for both reduced 
theories is that the ratio e of the gyroperiod to the macroscopic time scale is 
small. In the drift kinetic theory e is the same as the ratio of the gyroradius 
to macroscopic scale length while in the gyrokinetic theory small spatial 
variations are permitted. Because of the reduction of the phase space some 
information of the particle motion is missed. In this respect it may be noted 
that the reduced-phase-space kinetic theories are developed via expansions 
in e, the convergence of which is not guaranteed. This gives rise to the 
question: is the missing information important? To address this question we 
compared axisymmetric Vlasov equilibria near the magnetic axis with drift 
kinetic and gyrokinetic ones [25]. It turned out that because of missing the 
above constant of motion, C, in the latter case the on axis current density can 
differ appreciably from the "actual" Vlasov current density. Also, unlikely 
Vlasov theory, the reduced kinetic theories can not distinguish a straight 
from a curved magnetic axis. 

The aim of the present contribution is twofold: first to extend the above 
local constant of motion, C, to the more general class of helical symmetric 
equilibria and second to compare the Vlasov helical symmetric equilibrium 
characteristics near the magnetic axis with respective drift kinetic and gy- 
rokinetic characteristics. The new local constant of motion is derived in Sec. 
II. Then, Vlasov theory is compared near axis with reduced kinetic theories 
in Sec. III. Sec. IV summarizes the conclusions. 

II. A third Vlasov constant of motion near magnetic 

axis 

We will employ the following form of Vlasov equation [26] : 

| + vV/ + e r (E + vxB)^ + (e r vxVxv)^=0. (1) 

at ovi ovi 

To derive ([T]) one uses general orthogonal coordinates (xi,X2,Xs) with unit 
basis vectors = Vxj/|Vxj| (i = 1,2,3), expresses the "microscopic fluid" 
velocity in the basis of e, as v = and uses the "microscopic fluid" mo- 
mentum equation, 

<9v _ „ 

— + v ■ Vv = E + v x B, 

at 
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where E and B are the electric and magnetic fields consistent with Maxwell 
equations. The term "microscopic fluid" relates to the fact that the Vlasov 
equation is an approximation to the N-particles Liouville equation, which 
replaces the N particles by a continuum. This is sometimes termed "fluid 
approximation". To avoid confusion, however, the usual term "particle" 
will be adopted in place of "microscopic fluid" . Also, for the sake of notation 
simplicity and without loss of generality we will consider only ions and employ 
convenient units by setting m = q = c = 1 where m and q are the ion mass 
and charge and c is the velocity of light. 

In connection with the helical symmetry we specify the coordinates as 
(xi = r, x 2 = a<j)+(3z, £3 = — (3<f)+az), where a and (3 are parameters. Helical 
symmetry means that any quantity does not depend on x 3 . Translational 
symmetry and axisymmetry are then recovered as particular cases for (a = 1 , 
(3 = 0) and (a = 0, (3 = 1). The respective unit vectors are 

aes + (3re z are z - Pej, 

e x = e r , e 2 = — e 3 - 



[a 2 



/3 2 r 2 ) 1/2 ' (/3 2 + a 2 r 2 ) 1/2 ' 



Note that these basis vectors are in general non orthogonal. A helical mag- 
netic axis however is located at a constant distance r from the z-axis. Since 
we will study the equilibrium near the magnetic axis we set r = 1 thereon 
and make the choice a 2 + (3 2 = 1. Consequently, the above basis vectors re- 
duce to (e! = e r , e 2 = ae^ + (3e z , e 3 = —f3e ( p + ae z ) and become orthogonal 
on axis, with e 2 along and e 3 perpendicular to the axis. Thus, (JT|) can be 
applied in the neighbourhood of axis. The magnetic field near axis can be 
written in the form 

T3 I u Jae 2 -/3e 3 (3e 2 + ae 3 
B = -e, + B e z = - + B . 

By the choice a£>o = (31 we take on axis B = J/ae 2 as it should be. Also, 
we examine the case of E = on axis but V/ 7^ thereon. The latter 
assumption may be regarded as extraordinary because V/ is related to the 
density gradient on axis which for usual peaked density profiles is expected 
to vanish. The reason for assuming V/ 7^ on axis is that for axisymmetric 
equilibria, if ([1]) is solved near axis by the method of characteristics under 
the assumption V/ = on axis the usual toroidal momentum constant of 
motion, p^, is missed, while is recovered if V/ 7^ thereon |25j . 
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Using the above basis for the helically symmetric equilibria under consid- 
eration Eq. ([T]) near axis becomes 

df df df df df 

vi- \-v 2 - hwi- hfc \-w 3 — = 0, (2) 

OXi OX 2 OVi ov 2 ov 3 

where 

wi = (av 2 - f3v 3 ) 2 -v 3 -, w 2 = av x (/3v 3 - av 2 ) , w 3 = Vi~+Vif3 (av 2 - f3v 3 ) . 

a a 

(3) 

The characteristics of ([2]) are given by the solutions of 

dx\ dx 2 dv\ dv 2 dv 3 
dvi dv 2 dwi dw 2 dw 3 

Integration of the last equation yields the new constant of motion: 

d{(3 2 - a 2 ) (f3v 3 - av 2 ) +alhx\ [a 4 v 2 - a 2 /3 2 v 2 + a^ 3 v 3 -/?(/ + a 3 v 3 )] \ 

^ — v 3 H o 

(a 2 - f3 2 ) 2 

for a (5) 
n o (£{v 2 -v 3 f 

C = 2v 3 -v 2 -— for a = /3. (6) 

The respective axisymmetric constant of motion, C = v z + /In \v$\, is re- 
covered from (jSJ) for (a = = 0) and the translational symmetric con- 
stant of motion, C = v z , for (a = 0,/3 = 1). Also, conservation of energy 
H = 1 /2 (v 2 + v 2 + v 2 ) follows from the manipulation of equations 01]) pre- 
sented in Appendix. Distribution functions of the form f(H, C) could create 
a helical current on axis because of the dependence of C on v 2 additionally 
to that in the logarithmic term. This is different from the axisymmetric case 
in which toroidal axisymmetric currents on axis are not possible at all for 
the same class of distribution functions because of the dependence of C on 
\v^,\. This difference may be expected because for translational symmetric 
equilibria which are recoverable from the helically symmetric ones the respec- 
tive distribution functions f(H,C) = f(H,v z ) can create currents on axis. 
However, such a helical equilibrium is not physically acceptable because of 
the dependence of C on v 3 (but not on v\) implying the creation of a flow 
perpendicular to the axis. Similar unphysical flows are created in the ax- 
isymmetric case because of the respective dependence of C on v z (but not on 
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v r = v\). As discussed in jH] possible reasons of this unphysical behaviour 
are: (i) the lack of a potential fourth constant of motion involving V\ in 
such a way that creation of regular poloidal flows be possible, (ii) non taking 
into account here the MHD property of coincidence of the magnetic surfaces 
with the flow surfaces, (iii) additional drifts because of the curvature and 
torsion of the magnetic field which are eliminated in translational symmetric 
geometry and (iv) potential damping of unphysical flows in the framework 
of a collisional kinetic theory. The above consideration shows that straight, 
circular and helical magnetic axes are well distinguished in the framework of 
Vlasov theory. 

Unlikely the axisymmetric case [25] the generalized momentum constant 
of motion, p x3 = v 3 + A 3 , where A(xi,X2) is the vector potential, can not 
follow from @ though we have assumed V/ ^ on axis. This is obtained 
form the Lagrangian 



even in the presence of an electric filed on axis associated with the electro- 
static potential $. Also, a generic non local derivation off p x3 is provided in 



III. Comparison of Vlasov with drift kinetic and 
gyrokinetic theories. 

The near axis consideration of helically symmetric equilibria will be re- 
peated in the framework of the drift kinetic and gyrokinetic theories which 
were employed in Ref. [25] on an individual basis below. For convenience a 
brief review of these theories will also be given here because otherwise fre- 
quent reference to [25] would make reading of the present section tedious. 

Drift kinetic theory 

The drift kinetic theory established in Ref. [31] is based on the Little- 
john's Lagrangian for the guiding center motion [32] extended to include the 
polarization drift in such a way that local conservation of energy is guaran- 
teed. Introducing the modified potentials, A* and $*, this Lagrangian can 
be written in the concise form [33] 



L = o (^i + +i 2 + &l) + x 2 ) + x ■ A(xi,x 2 ) 



[27\ 




L = A* • R — $*, 



(7) 
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where 



A* = A + v\\b + V E , (8) 
= $ + pB+±(tf\+v%), (9) 

v, = ^ do) 

Here, $ and A are the usual electromagnetic scalar and vector potentials 
and b = B/.B. The drift kinetic equation for the guiding center distribution 
function f(R.,vn,fJ,,t) (with fi — 0) acquires the form 

^-+V-V/ + ^ = 0. (11) 

OT OV\\ 

By means of the Euler-Lagrange equations following from (j7j) the guiding 
center velocity, V, and the "acceleration" parallel to the magnetic field, v\\, 
can be expressed as: 

v + (12) 

B \\ B \\ 

"17* T3* 1 

where 

E * = — ~ m- B " = VxA *- (14) 

Sj] r = B*.b = S + «||b- Vxb + b- VxV E . (15) 

Explicit expressions for V and V\\ are given by Eqs. (3.24) and (4.17) of Ref. 
[3T] . Also, it is noted here that Eqs. ( lT2"j) and f[T3^) have similar structure as 
the respective gyrokinetic equations of Refs. [38] and [39] [Eqs. (5.39) and 
(5.41) therein]. Since B*\ appears in the denominators of (TT2T) and (jTBl a sin- 
gularity occurs for B* = 0. For E = this singularity can be expressed by the 
critical parallel velocity v c = — fi/(b ■ Vxb), where Q is the gyrofrequency. 
Therefore, the theory is singular for large \v\\ \ at which V and v» diverge and 
consequently non-casual guiding center orbits occur and the guiding center 
conservation in phase space is violated [30]. It is the vy-dependence of A* 
[Eq. OH])] that produces the singularity. In order to regularize the singularity 
v\\ in OH]) can be replaced by an antisymmetric function g(z) with z = vu/vq, 
where vq is some constant velocity [30] [3T] |33] . The nonregularized theory 
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presented here for simplicity is obtained for g{z) = z. In the regularized 
theory g(z) ps z should still hold for small \z\. For large \z\, however, g must 
stay finite such that with v ^> ^thermal one ^ as v o9(°°) < v c . A possible 
choice for g(z) is g(z) = tanhz. 

Since for a helically symmetric equilibrium the magnetic field on axis 
becomes purely helical and dependent only on r, one readily calculates near 
axis with E = thereon: 



Vxb = Vxe 2 = -e 2 , B* = B + v\\-e z , ffi = B + vn-P (16) 

/y H /y II 'I ry> 



V5(r) = i, & = pB + ±v?>, E* = -^e r (17) 
dr 2 11 ar 

and consequently 

B v\a n dB 

v = l W 2 + 7b**' - -5* IF e " (18) 

t)|| = 0. (19) 

As expected on axis the guiding center velocity consists of a component 
parallel to B and the curvature and grad-5 drifts. Also, the "acceleration" 
U|| vanishes because there is no parallel force and the drift kinetic equation 
(fTTT) becomes 

v 'lh - < 20 > 

Because of (Tl9l) vu is a constant of motion and therefore distribution functions 
of the form f(H,v\\), where H = [iB + l/2t|, either can or can't produce 
helical currents by choosing / either odd or even function of v\\. This prop- 
erty holds also for axisymmetric and translational symmetric equilibria fZ5\ 
because f fl9|) keeps valid irrespective of the kind of symmetry. Thus, un- 
like Vlasov theory the drift kinetic theory can not distinguish equilibria of 
straight, circular or helical magnetic axes. Also, according to the results of 
section III the respective obtainable near axis Vlasov class of distribution 
functions f(H,C) can not produce physically acceptable currents on axis^. 
Consequently, since near axis the overwhelming majority of the particles are 



Physically acceptable Vlasov currents can be produced by the class f(H,p x3 ) which, 
however, corresponds to drift kinetic distribution functions of the form f(H,p^. 3 ), where 
oi 3 = A 3 +v r 



near axis p x3 = A3 + W11&3, as it follows from the Lagrangian (0. 



s 



passing, the parallel currents constructed in the framework of the drift ki- 
netic theory may differ from the "actual" ones. Also, unlike in the case of 
axisymmetric equilibria [25] , the B*\ -singularity is present as indicated by Eq. 
(US]). 

Gyrokinetic theory 

We will use the gyrokinetic equations of Ref. [35] which have been em- 
ployed to a variety of applications (see for example Refs. [IQJ HTJ 02]). Eq. 
(jTTj) remains identical in form where /(R, vu, fi, t) is now the gyrocenter dis- 
tribution function for ions. The gyrocenter velocity and "acceleration" are 
given by 

v = T,|,b + |r (v e + v Vb + V c ) , (21) 

o|| 

v\\ = -— V • (v¥ + fiVBo) . (22) 

Here, B is the equilibrium magnetic field, b = B Q /B , 

B* l{ = (B + <U||Vxb ) -b , 

$ stands for the perturbed gyroaveraged electrostatic potential, and the 
ExB-drift velocity Ve, the grad--B drift velocity V^ B , and the curvature 
drift velocity V c are given by 

Ve _v™ (23) 

Vvs-^BoXVBo, (24) 

V -^| b » XV ( P » + f)' < 25) 

Note that as in the case of drift kinetic theory a similar singularity occurs 
at Bq» = 0. In numerical applications this singularity was "avoided" by 
approximating Bq» = B (see for example Refs. [IQIIIT]). Consideration of 
the above equations for a helically symmetric equilibrium with E = near 
axis yields relations similar to ( I18j) . (TL9l) and ( 1201) . Therefore, the above found 
discrepancies of the drift kinetic theory with the Vlasov one persist in the 
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framework of the gyrokinetic theory. The structure of the reduced kinetic 
equations in conjunction with the symmetry of the equilibrium considered 
clearly indicate that this conclusion is independent of the particular drift 
kinetic or gyrokinetic equations. 

IV. Summary 

We have found a new constant of motion near the magnetic axis of heli- 
cally symmetric equilibria in the framework of Vlasov theory [Eqs. (j5J) and 
flSJ)]. On account of this constant of motion a comparison of the Vlasov 
equation with either the drift kinetic or the gyrokinetic equation indicates 
that the current densities near the magnetic axis in the former case may be 
different from the ones in the latter case. Also, unlike Vlasov, the reduced 
kinetic theories can not distinguish equilibria with straight, circular or helical 
magnetic axes. This discrepancy is due to the loss of new Vlasov constant of 
motion in the reduced phase space thus indicating that this reduction, even 
if made rigorously so that local conservation laws and Liouvillean invariance 
of the volume element is guarantied, is associated with the loss of nontrivial 
physics. This could put certain limits on the conclusions from drift kinetic 
or gyrokinetic simulations for time scales large compared to the gyroperiod 
which is the case of an equilibrium. Though the derivations of both theories 
are formally correct the convergence of the e-expansions is presumably not 
uniform for all times which could lead to the noticed discrepancies calculated 
in this study. 

In addition, a singularity which occurs in both drift kinetic and gyroki- 
netic theories for large parallel particle velocities, present for helically sym- 
metric equilibria, is usually eliminated in the literature by a rough approxi- 
mation. Alternative to the rather artificially imposed regularization reported 
in Sec. Ill, one might remove this singularity in the integrals associated with 
moments of the drift kinetic equations (see the generic relations (8.14)-(8.17) 
of [21] for the self consistent charge, current, energy and energy flux densi- 
ties) by Cauchy principle value integration provided that such an integration 
would be physically justifiable. This requires further investigation. 

Appendix: Energy conservation near axis 

From 

dx\ dv\ dx 2 dv 2 
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[see Eqs. (j3J] it follows 



Also, employing 



one has 



Widxi + w 2 dx 2 = -d(v{ + v 2 ) (Al) 
2 

dx\ dx2 dv 3 
Vi v 2 w 3 

w 1 dx 1 + w 2 dx 2 = (wxVi + 1^2) — (A2) 

w 3 

Explicit calculation by means (J3J) yields 

WiVi + w 2 v 2 = -w 3 v 3 . (A3) 
From (TAT]), (jM]) and (TA3]) it follows 

-{v\ + w 2 + U3) = const. 
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